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Abstract 

The propagation of an electrostatic wave packet inside a collisionless and initially Maxwellian 
plasma is always dissipative because of the irreversible acceleration of the electrons by the wave. 
Then, in the linear regime, the wave packet is Landau damped, so that in the reference frame moving 
at the group velocity, the wave amplitude decays exponentially with time. In the nonlinear regime, 
once phase mixing has occurred and when the electron motion is nearly adiabatic, the damping 
rate is strongly reduced compared to the Landau one, so that the wave amplitude remains nearly 
constant along the characteristics. Yet, we show here that the electrons are still globally accelerated 
by the wave packet, and, in one dimension, this leads to a non local amplitude dependence of the 
group velocity. As a result, a freely propagating wave packet would shrink, and, therefore, so would 
its total energy. In more than one dimension, not only does the magnitude of the group velocity 
nonlinearly vary, but also its direction. In the weakly nonlinear regime, when the collisionless 
damping rate is still significant compared to its linear value, this leads to an effective defocussing 
effect which we quantify, and which we compare to the self-focussing induced by wave front bowing. 
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I. INTRODUCTION 



Collisionless dissipation was first described by Landau in his famous paper Ref. jl| pre- 
dicting tlie damping of an electrostatic wave, and is due to tlie global acceleration of the 
charged particles by the electrostatic potential, as was later discussed in Refs. . Landau 
damping is certainly one of the most important and basic features of wave propagation in 
a plasma, it has many applications (e.g. current drive in a tokamak, j4, O]), and has been 
extensively discussed in the scientific literature, which eventually led to its rigourous proof 
(see Ref. jol and references therein). It is nevertheless not ubiquitous since, as shown by 



O'Neil in Ref. an electron plasma wave (EPW) with a large enough initial amplitude 
remains essentially undamped, because the rate of energy exchange between the electrons 



and the EPW decreases on a timescale of the order of the bounce period, Te = 271 ^Jm/ eE^k, 
m being the electron mass, — e its charge, Eq the EPW amplitude and k its wave number. 
Then, if lob = ^t^/Tb is much larger than the Landau damping rate, vl, the EPW amplitude 
remains nearly constant (see Ref. O'Neil's theory was the subject of some controversies 
(see Ref. js] and references therein) but was checked experimentally in Ref. js] where it was 
shown that, if an EPW was excited with a small initial amplitude, ub < i^l, it would damp 
at the rate predicted by Landau while, if ub ^ i^l, after a few oscillations at a period close 
to Tb, the EPW amplitude would only slowly decrease on a timescale much larger than uj^^ 
due to some finite dissipation from the detection equipment. 

The nonlinear reduction of collisionless damping was recently found to have important 
implications in inertial confinement fusion, using the indirect drive approach 10|, since it 
could lead to Raman reflectivities much larger than linear theory would predict on the 
grounds of Landau's prediction for the EPW damping rate. This was unambiguously shown 
experimentally in Ref. and these experimental results were qualitatively reproduced 
using kinetic simulations by several authors (for example in Refs. 12|, ll3|). The nonlinear 



reduction of the collisionless damping rate of an SRS-driven plasma wave was furthermore 



derived theoretically in Refs. [IJ, ll5|; two theories supported by Vlasov simulations and 



essential 
in Refs. 



y equivalent over a finite range of wave amplitudes, as shown in Ref. [16| . Moreover, 
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by comparing results found in Vlasov simulations to those derived using a 
nonlinear envelope code following the theory of Ref. {iqI, levels of Raman reflectivity way 
beyond those predicted by linear theory were very clearly associated with the vanishing of 
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the collisionless damping rate. Finally, anomalously large reflectivity levels measured at the 



10| . and analyzed in Ref. 20|, seem to be due to nonlinear kinetic 



National Ignition Facility 
effects. 

Now, as a driven EPW grows in an initially Maxwellian plasma, it keeps on increasing 
the electron kinetic energy in an irreversible fashion, because it keeps on trapping more and 
more electrons. Moreover, the number of trapped electrons quickly increases with the EPW 
amplitude, and so does the amount of electrostatic energy which is dissipated. This seems 
to contradict the previous discussion on the nonlinear reduction of the collisionless damping 
rate which should vanish after a few bounce periods, and therefore more rapidly if the wave 
amplitude is large. One of the scopes of this paper is to resolve this apparent paradox by 
showing that an undamped wave may still be subjected to collisionless dissipation, which 
manifests itself in the nonlinear variations of the EPW group velocity, Vg. In one dimension, 
Vg decreases in that space region where the wave amplitude increases along its direction 
of propagation, and remains almost constant elsewhere. It is quite clear that this would 
automatically reduce the total energy of a freely propagating and undamped wave packet, 
simply because the total area located under the wave envelope would shrink, until a shock- 
like situation would occur (see Fig. [2]). Understanding this point allows us to discuss 
physically, in Section |Tll the results previously found in Ref. 2l|] on the nonlinear group 
velocity of an EPW. In particular, we clearly specify when, and how, collisionless dissipation 
enters to let us derive an envelope equation that does not follow from the variational approach 
developped by Whitham in Ref. 22|. Moreover, one important point of this paper is to 
show how the reasoning made in one dimension straightforwardly generalizes to allow for 
a mult i- dimensional geometry. This lets us predict that collisionless dissipation leads to a 
non zero transverse component of the EPW group velocity that was quite unexpected, and 
which we very precisely quantify. In particular, we discuss its ability to offset, or to locally 
reinforce, the focussing effect due to the wavefront bowing induced by the EPW nonlinear 
frequency shift (see for example Ref. 13|) for a plasma wave resulting from stimulated 
Raman scattering. 

This paper is organized as follows. The next Section is devoted to one-dimensional 
wave propagation, starting with the situation when the wave amplitude is uniform. In 



this situation, we physically discuss the envelope equation derived previously in Refs. 
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19| . and emphasize the very particular role played by the trapped electrons, which leads 
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to an equation that would not follow from Whitham's theory. This envelope equation is 
then straightforwardly generalized to allow for a inhomogeneous wave amplitude, which 
lets us derive in an extremely simple fashion the results of Ref. 21| on the EPW group 
velocity. In particular, it is unambiguously shown that the nonlinear variations of Vg allow 
for a collisionless dissipation mechanism of the electrostatic energy different from Landau 
damping. In Section llllt the results obtained in one dimension are very easily generalized to 
multi-dimensional wave propagation, leading to new nonlinear focussing/defo cussing effects 
which were, in some way, qualitatively discussed in Ref. 23j, and which we quantify. This 
allows us to make comparisons with the focussing effect due to wave front bowing, and to 
discuss which of these effects prevails, depending on time and on space location. Section [IV] 
concludes and summarizes our work. 



II. ONE-DIMENSIONAL WAVE PROPAGATION 
A. Uniform wave amplitude 

In this Subsection, we consider a nearly monochromatic electron plasma wave, with fre- 
quency u and wave number k, whose amplitude varies slowly in time while remaining space 
independent. We moreover restrict here to the situation when the time it takes for an elec- 
tron to complete one "frozen" orbit (i.e., corresponding to a fixed wave amplitude) is much 
smaller than the typical timescale of variation of the EPW amplitude. For a harmonic wave, 
this condition is achieved when ub ^ 7, where 7 is the EPW growth rate. In this case, 
one may assume that the electron orbits are infinitely close to the "frozen" ones, and that 
these orbits are completed in an infinitely small time (which would correspond to making 
use of the adiabatic approximation). In other words, the electron motion is "enslaved" to 
the variations of the wave amplitude. Then, in this situation, it is clear that any observable 
related to the electrons (energy, momentum, charge density,. . . ) would vary at a rate nearly 
proportional the the wave growth rate, 7. As for the ions, for the sake of simplicity, we 
henceforth assume that they are immobile. 

Let us now consider an EPW whose amplitude, Ep, grows from -Emin to -Emax, and then 
decreases back to -Emin- If the electron motion has kept on being adiabatic during these 
variations of Ep, it is easily shown that the total energy of the electrons which have never 
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been trapped is only a function of the EPW amplitude. Hence, when the wave grows, the 
gain in energy of the untrapped electrons is reversible and given back to the wave as it decays. 
Now, whatever Av, adiabatic electrons with initial velocities f(^± Af , where v^ = uj/k is the 
EPW phase velocity, are all trapped at the same time, and their trapping leads to a jump 
in the kinetic energy by a quantity proportional to 



A£ = m[fo{y^ - Av) - fo{v^ + Av)]v^Av, 



(1) 



where /o is the electron distribution function in the limit of a vanishing wave amplitude. 
Note that AS is non zero only because fo{v^ — Av) 7^ /o(i'(/> + At>) and, actually, AS is 
positive when /q is a decrea sing function of velocity, as is the case for a Maxwellian. Now, as 
shown theoretically in Ref. [l9|, and numerically in Refs. 2J,|25|, while Ep is decreasing, the 
electrons are detrapped nearly symmetrically with respect to the phase velocity (see Fig. [1] 
for a sketch of symmetric detrapping). Consequently, detrapping would lead to a change 
in the electron kinetic energy by AS' = —AS but, now, with fo{v^ — Av) = f^iy^ + Av). 
Hence, AS' = 0, the energy gained by the electrons through trapping is not given back to 
the wave when they are detrapped. We therefore conclude that, when the electron motion 
is nearly adiabatic, only trapping may lead to an irreversible increase of the kinetic energy, 
and therefore to the collisionless dissipation of the electrostatic energy. This point is actually 
not new and was already discussed in Refs. [23|, l26|. However, the scope of this paper is to 



describe in detail the implications of collisionless dissipation in the limit when the electron 
motion is nearly adiabatic, and to make clear how it may differ from Landau damping, which 
is one distinction we could not find in any of the papers we know of. Our discussion will be 
based on the envelope equation for the EPW amplitude, Ep. 



Such an equation was derived by Yampolsky and Fisch in Ref. 15| for a growing driven 
wave, and these authors found 



dujXr{dt + v)Ep = Ed, 



(2) 



where Xr is the real part of the electron susceptibility, E^, accounts for the effect of the 
drive, while u is directly related to the rate of energy gained by the electrons. As discussed 
before, for adiabatic electrons, the rate of energy variation is proportional to the EPW 
growth rate, and therefore so should be u in the limit when ^ 7. This is exactly the 
result found by Yampolsky and Fisch in Ref. [l5| so that, when ub ^ 7, one may write 
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FIG. 1: (Color online) Space portrait of electrons with initial distribution function fi = Ni6{v — 
V(fi — At>) + N26{v — v^ + Aw) acted upon by an electrostatic wave, (a) The wave amplitude is zero, 
(b) The wave has grown slowly enough to induce adiabatic motion and has trapped the electrons 
which, due to action conservation, all lie on the same orbit, (c) In the limit when the period of 
the trapped orbit is infinitely small compared to the typical timescale of variation on the EPW 
amplitude, when this wave decays, half electrons are detrapped with v > and half with v < v^. 
Then, again due to action conservation, when the EPW amplitude has decreased back to zero the 
distribution function is /o = [(A''i + N2)/2][6{v — — At>) + 5{v — V(f) + At>)]. 

udujXrEp = duiX'r^tEp, where Xr is a dimensionless function of the EPW amplitude, and 
Eq. ([2]) becomes 

[d^Xr + d^Xr)dtEp = Ed. (3) 



Unlike in Ref. 15|, in this paper we make the choice to no longer term u the "nonlinear 
Landau damping rate of the EPW" when Eq. ([2]) is essentially equivalent to Eq. ([3]). In- 
deed, when this condition is met and when the EPW is freely propagating [E^ = 0), its 
amplitude remains constant, the wave is undamped. Yet, if the wave is driven and keeps 
on trapping new electrons, it increases the electron energy in an irreversible fashion, and is 
therefore subjected to collisionless dissipation, which is accounted for by the term d^jXr^tEp 
in Eq. (|3]) . This is a first illustration of the distinction we make between Landau damping 
and collisionless dissipation. 



It is very interesting to note t 
obtained by Benisti et al. in Refs. 



lat the same kind of envelope equation as Eq. ([3]) was 
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19| without ever calculating the rate of energy gained 
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by the electrons. Indeed, starting from Gauss law, Benisti et al. straightforwardly find 
XiEp = Ed, where Xi is the imaginary part of the electron susceptibility. Moreover, using an 
argument resting on the orbit symmetry, these authors show that, in the limit cj^ ^ 7, the 
contribution to Xi of ths trapped electrons is negligible, so that Xi ~ xT^^^^ where Xi^^'^ is the 
value of Xi obtained by accounting only for the untrapped electrons. For a slowly varying 
wave amplitude, it is easily shown that Xi^^'^ is nearly proportional to the wave growth rate, 
which leads to the following envelope equation, 

d^xT^'dtEp = Ed, (4) 

where Xr^^^ is the sole contribution of the untrapped electrons to the real part of the suscepti- 



bility. Since, as shown in Ref. 16|, the results obtained by Benisti et al. and Yampolsky and 
Fisch match over a range of wave amplitudes where the condition Ub ^ 7 holds, Eq. (j4]) may 
be identified with Eqs. ([2]) and ([3]). This shows that, when ub ^ 7, i^/'j ~ —duiXl^/duiXr, 
where xl^ is that part of Xr only due to the trapped electrons. This formally relates colli- 
sionless dissipation to trapping, as expected from our previous discussion. 

The physics content of Eq. (|2]) seems extremely clear, the rate of variation of the wave 
amplitude results from the balance between the energy gained from the drive and that given 
to the electrons. When this equation may be cast in the form of Eq. (j3]) or Eq. (|4]) then. 



as shown in Refs. 
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21| . duiXr^^^ is larger than d^iXr, all the more as the (linear) Landau 
damping rate of the EPW is large. This just reflects the fact that a driven EPW would grow 
slower if a large fraction of the energy it gains from the drive is given back to the electrons 
through trapping. 

Note, though, that the theory by Yampolsky and Fisch [15], based on energy conservation, 
only holds when the EPW keeps growing. By contrast, the argument on the symmetry of the 



trapped orbits given by Benisti et al. in Ref. 19|] is clearly valid whether the wave grows or 



decays. Hence, provided that the electron motion is nearly adiabatic, Xi ~ Xf^^^ ~ id^^Xr^^^^ 
so that Eq. @ holds independently of the sign of dtEp. Nevertheless, although this equation 
is always formally valid, for a given Ep < -Emax, d^Xr^^^ does not assume the same values 
when Ep has kept on growing as when Ep has previously decayed from -Emax- Indeed, as 



shown in Ref. 



19| . due to symmetric detrapping, the contribution to d^^Xr^*^ of the detrapped 



electrons is negligible. Hence, Eq. (j4j) would be better written 

d^xfdtEp = Ed, (5) 



where xf^ is calculated by accounting only for the electrons which have never been trapped. 
The envelope equation for the plasma wave is therefore non-local in the EPW amplitude since 
the value of d^xt^ at time t does not only depend on Ep(t) but also on maxt/^t[Ep(t')]. As 
shall be seen in the next Subsection, this non-locality entails the reduction of the total energy 
of a freely propagating wave packet. Moreover, when the wave is essentially undamped (i.e., 
when its amplitude remains constant along the characteristics), this is the main cause for 
the dissipation of the electrostatic energy. 

From the previous discussion we conclude that, if one were to write a Lagrangian for the 
wave-particle interaction, this Lagrangian would be non-local in Ep, so that it would not be 
quite easy to derive the envelope equation for the EPW from a variational approach, and 
this program will not be pursued in this paper. However, it is very interesting to discuss 
physically the correction that needs to be made to the variational approach developed by 
Whitham in Ref. 22| in order to recover Eq. ([5]). In Whitham's theory, the Lagrangian is 
L{Ep) = J^^ll + Xr{E')]E'dE' (notice that this is the Lagrangian one would derive within 
Whitham's formalism which does not account for dissipation. This is not the Lagrangian 
that should be used for the nonlinear propagation of an EPW we are addressing here). Then, 
for a freely propagating wave whose amplitude is uniform, since Whitam's theory cannot 
account for dissipation, it would trivially predict d^jXrdtEp = 0. For a driven wave, this 
equation would generalize into d^jXrdtEp = Ed. Based on this, we now rewrite Eq. ([5]) as 

d^XrdtEp + d^lxf - xAdtEp = E^. (6) 

The second term in the left-hand side of the previous equation clearly is the correction 
that must be made to Whitham's theory in order to accurately describe the nonlinear time 
evolution of the EPW amplitude, and its physics content is transparent. Indeed, {xf^ ~ 
Xr) is the contribution to the real part of the susceptibility of those electrons which have 
been trapped at least once, while dtEp is proportional to the trapping rate. Hence, the 
correcting term is directly related to trapping, because in the limit ub ^ 1 only trapping 
entails collisionless dissipation. Understanding this will let us generalize Eq. ([6]), in a very 
simple and powerful fashion, in order to allow for a multi-dimensional variation of the wave 
amplitude, as we shall show it in the remainder of this paper. 

Before this, we would like to stress that collisionless dissipation, and not trapping by 
itself, is at the origin of the correcting term in Eq. ([6]). Indeed, consider again an EPW 
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whose amplitude grows from -E'min to -E'max, and then decreases back to i?min- If the initial 
distribution function is symmetric with respect to on a given velocity interval, 1^,, then, 
from Eq. ([T]), the trapping of electrons with initial velocities in does not change the total 
kinetic energy. Moreover, those electrons, either trapped or not, do not contribute to d^^^Xr- 
Hence, as long as only those electrons have been trapped, there has been no irreversible 
increase in the electron kinetic energy, and the correcting term in Eq. (|6]) is strictly zero 
although trapping did occur. 



B. Non-uniform wave amplitude 



In the previous Subsection, we concluded that Whitham's result had to to be corrected 
to account for coUisionless dissipation which, when the electron motion is nearly adiabatic, 
is only due to trapping. Using this point, it is extremely easy to generalize Eq. (E]) to a 
wave whose amplitude only varies along its direction of propagation, x. Indeed, to do so, 
we only need to add to the left-hand side of Eq. (EJ the term —dkXrdxEp coming from 
Whitham's theory 2^], and to express the fact that the correction to Whitham's result is 
proportional to the rate of electron trapping, and therefore proportional to the EPW growth 
rate calculated in the frame moving at the phase velocity with respect to the laboratory frame. 
Hence, this term is d^[xf ~Xr][dtE, + v^d,Ep] so that, once Landau damping has vanished, 
the nonlinear envelope equation of the EPW is 



dxr dEp dxr dEp ^ d[x 



cff 
r 



Xr 



dEp 
dt 



dEp 
dx 



E.. 



(7) 



dtu dt dk dx doj 

This is exactly the result found in Ref . {21!] , which was derived in the Appendix of that paper 
directly from Gauss law and not by using a variational approach. In Ref. 2l|, we already 
discussed physically that our result differed from the one derived by Whitham because of 
coUisionless dissipation. Here, we will push the reasoning further and we will show in a few 
lines that the nonlinear variation of the EPW group velocity, as derived in Ref. 21|, does 
entail the coUisionless dissipation of the electrostatic energy. Here, we call 

_ dkXr , d^xf-Xr] 



,cff 



(8) 



the EPW group velocity because, if the wave was freely propagating, its amplitude would 
remain constant along the lines x = Vgt From the results of the previous Subsection, we 
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FIG. 2: (Color online) Space profile of the plasma wave (in arbitrary units), at time ti (black solid 
line) and at time t2 > ti (blue dashed line), when the group velocity of the wave packet (whose 
amplitude is indicated by the arrows) decreases with the EPW amplitude at the rear side, and 
remains fixed at its minimum nonlinear value at the front side. One clearly sees that the wave 
packet at time t2 is narrower than at time ti. 

clearly understand that Vg differs from —dkXr/dujXr because of the collisionless dissipation 
entailed by trapping. 

Now, in order to discuss the implications of the nonlinear variations of Vg in the decrease 
of the total electrostatic energy, we first consider a freely propagating wave packet whose 
evolution is given by Eq. ([7]) with = (we therefore consider here a strongly nonlinear 
regime, where we neglect the effect of damping which will be discussed below, and which 
is assumed to be localized in a space region where the EPW amplitude is very small com- 
pared to its maximum value). Since the wave packet is essentiall y u ndamped, its maximum 



amplitude remains nearly constant. Moreover, as shown in Ref. [2l|, Vg = — 2/(kd^Xr)i 



so that Vg < v^. Hence, along the rear side of the wave packet, i.e. before its maximum. 



Ep{x + Vfft) keeps on increasing so that, as shown in Ref. |2l|, Vg decreases (for the sake of 
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simplicity, we assume here that the wave packet has no local maxima). In front side of the 
wave packet Ep{x + v^t) keeps on decreasing, so that r[iax.ti <t{Ep[x — v^{t — t'),t']\ remains 
constant and so does the number of electrons which have been trapped by the EPW. Since 
d^X^^ essentially depends on these parameters, the group velocity remains nearly fixed at 
its minimum nonlinear value along the whole front side of the wave packet. Hence, the rear 
side moves faster than the front side so that the wave packet shrinks while its maximum 
amplitude remains constant, as illustrated in Fig. [2l which automatically reduces the total 
electrostatic energy, W = E^dx. This goes on, at most, until a shock would occur in 
the wave profile. At this point, the EPW amplitude, as seen by the electrons, would no 
longer vary slowly so that the adiabatic approximation would break down and our analysis 
would no longer hold. What happens then is outside the scope of this paper and will not be 
addressed here. 

Our prediction on the shrinking of the wave packet is consistent with the numerical results 



of Ref. [26|. However, that paper only discusses the inhomogeneity of the Landau damping 
rate, larger at the rear side where the amplitude should therefore decrease more rapidly (or 
increase less rapidly for a driven wave) than at the front side, which leads to an effective 
shrinking of the wave packet. Here we show that another effect, the inhomogeneity of the 
EPW group velocity, which is one manifestation of collisionless dissipation but which occurs 
only after Landau damping has nearly vanished, may also reduce the space extent of the 
electrostatic wave packet. 

In order to clarify the role of each of these effects, let us consider a freely propagating 
EPW whose envelope equation is (see Ref. {l8|) 

dtEp + Vgd^Ep + vEp = 0, (9) 

where v is what we call the nonlinear Landau damping rate of the EPW. Using an integration 
by parts, one easily finds 

dt / E^dx = -2 vEldx + / d^VgE^dx. (10) 



As shown in Ref. [l^, t/ is a decreasing function ofY = J^^ujB{x — v^t' , t')dt' . Since > Vg, 
V is indeed larger at the rear side of the wave packet than at the front side as argued in 



Ref. 



. Actually, in Refs. we showed that i>{x,t) and Vg{x,t) remained very close 



to their linear values where Y{x,t) < 6, making, there. Landau damping the main cause 
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of collisionless dissipation. In the region where Y{x,t) ^ 6, v{x,t) quickly decreases with 
X, while Vg quickly increases towards v^. Hence, in the small space region where Y Q 
dissipation is, again, only due to Landau damping while the effect of d^Vg is anti-dissipative. 
Actually, in the weakly nonlinear regime when Y{x,t) ^ 6 close to the location of the 
maximum EPW amplitude, the variations of Vg would mostly entail the defocussing of the 
wave packet, in contrast with the strongly nonlinear regime illustrated in Fig. [2l Finally, 
wherever F > 6, z/ is strongly reduced compared to its linear value while, as discussed above, 
dxVg < in the rear side of the wave packet and, due to symmetric detrapping, dxVg in 
the front side. Hence, wherever V > 6, the nonlinear and non local variations of Vg are one 
main cause of dissipation, clearly more effective for longer and more intense wave packets. 

We therefore conclude that Landau damping, i.e. the decay of the EPW amplitude along 
the characteristics, and the decrease of the EPW group velocity along the direction of 
propagation of the wave packet, are two distinct manifestations of collisionless dissipation, 
which are effective on distinct space domains. Moreover, the variations of Vg may locally be 
anti-dissipative, which is most sensitive in the weakly nonlinear regime when Y Q close to 
the location of the EPW maximum amplitude. 

III. MULTI-DIMENSIONAL WAVE PROPAGATION 
A. Nonlinear envelope equation 

In this Paragraph, we use the same procedure as in the previous Section in order to 
generalize Eq. ([6]) to a multi-dimensional geometry. To do so, we just need to notice that, 
in Eq. ([6]), the term correcting Whitham's theory is proportional to the trapping rate, 
and, therefore, to the rate of variation of the EPW amplitude, calculated in the frame 
moving at the phase velocity with respect to the laboratory frame, and as seen by the 
electrons. When all electrons has the same transverse velocity, Vj_, this term is nothing but 
[dujX'f^ ~dijXr][dtEp + Vff,dx^^Ep + v_i_.'V _LEp], where, for paraxial wave propagation, x// is along 
the local direction of the wave vector, and is the gradient in the direction transverse to 
k. The correcting term is easily generalized to allow for a warm distribution of transverse 
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velocities so that, in more than one dimension, Eq. (Q becomes 



2 



where f{v±) is the distribution of transverse velocities. Notice that, even when f{v±) is 
isotropic, the third term in the left-hand side of Eq. (!TT|) may be non zero because d^Xr^ is a 
function of v±. Indeed, dujxf^ depends on the maximum wave amplitude, -Emax, experienced 
by the electrons, and -Emax is generally different, and usually larger, for the electrons exiting 
the wave packet than for those entering it. Hence, because the electron response to the wave 
is non local, the EPW group velocity may have a non zero transverse component, Vg^, which 
we now quantify in connection with the values assumed by the coUisionless damping rate, u. 



In Ref . 



14|, we showed that electrons such that Y3D = (^b{x// —v^t', x± — vj_t', t')dt' > 6 



i.e. those which have completed about one bounce period, no longer contribute to u, while 
all the others induce a nearly linear Landau damping of the EPW. This leads to the following 
estimate for u, 

where in the linear Landau damping rate and 'H(Y3d) is a Heaviside-like function, 
H(r3D) ^ if > 6 and ^(^30) ^ 1 if ^30 < 6. Similarly, in Ref. (l^ we showed 
that S^Xr*^ significantly departed from its linear value, d^^Xr^, only once Isd ^ 6, and that 
duXr^iXsD ^ 6) ^ di^Xr^- Moreover, when Y3D > 6, d^xt^ is essentially a decreasing func- 
tion of E^3,x{x//,xj_,t) = ma.xt> <:t[Ep{x // — v^t',x± — v±t',t')]. Using these results, we can 
now discuss the amplitude and orientation of Vg^^. 

In regions where the wave amplitude has been so weak that Y3D < 6 for all electrons, 
u ~ z/iin and d^ixf^ = c^wXr™ independently of v±, so that Vgj_ = if f{v±) is isotropic. 

When the maximum wave amplitude has reached a large enough value, there may be space 
locations where 1^313 > 6 for the electrons exiting the plasma wave packet while I3D < 6 for 
those entering it. There, for any Vj_ such that v±.'VEp < 0, d:^x'f^{^±) ^ d^jXr^i^vi.), 
and Vg^ oc / d^xf^ f {^i.)^±dv± is clearly oriented towards the outside of the wave packet, 
which would tend to make it self-defocus. Moreover, in such space locations, about half 
of the electrons (those exiting the wave packet) do not contribute to the Landau damping 
rate and, therefore, u ~ z/iin/2. Note that, as the wave grows, the region where Y3D < 6 for 
entering electrons is located further away from the center of the wave packet where the EPW 
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amplitude is largest, so that only the edge of the plasma pulse is affected by the defocussing 
induced by Vg^^, which, therefore, becomes less effective. 

In space domains where I3D > 6 for all electrons it is clear that u ^ 0. Moreover, since 
Ema.x is larger for electrons exiting the wave packet than for those entering it, and since, 
when > 6, d^xf decreases with E^,., d^xfiv^) < ^.xf (-t?^) if v^-^E, < 0. Hence, 
in such space regions, Vg^ is oriented towards the inside of the wave packet and would tend 
to make it self-focus. However, since d^^Xr^ varies slowly with -Emax, the magnitude of Vg^^ 
is usually weaker when it induces a focussing effect than when it tends to make the wave 
packet self-defocus. 

Let us now stress that Vg^^ is easily connected to collisionless dissipation since it is non 
zero just because, due to symmetric detrapping, di_jXr^ i^ a non local function of the EPW 
amplitude. As for collisionless dissipation, in the regime when 7, it is also only due 

to symmetric detrapping, as discussed in the previous Section. Actually, from Eq. (fTTj) one 
easily finds 

^ j E^r = -2 j vE^r + j d.^^Vg^^E^r + j V^-Vg^E^r, (13) 

where Vg^^ is the component of the group velocity along the x// direction [given by the sec- 
ond term in the lest-hand side of Eq. ([9])]. The effect of dx/^Vg^^ on the time variation of 
the electrostatic energy is exactly the same as in the one- dimensional situation. For each 
(x_L,t), wherever Y3Y){x//, x±,t) < 6, d^^^Vg^^ ^ 0, wherever Y^rtix//, x±,t) > 6, d^^^Vg^^ < 
in the rear side of the wave packet and dx^^Vg^^ in the front side, and in the small 
domain where Y3Y){x//,x±,t) ^ 6, dx//Vg^^ > 0. Hence, the effect of d^^^Vg^^ is mostly dissi- 
pative in the strongly nonlinear regime, for long and intense pulses, and anti-dissipative in 
the weakly nonlinear regime when Y3D ~ 6 close to the maximum wave amplitude. It is 
noteworthy that the same conclusion holds for each ^ ±-Vg^{v±), where Vg^{v±) is the con- 
tribution to Vg^ of those electrons which all have the same transverse velocity, v±. Indeed, 
for each {x//,t), wherever Y3r,{x//,x^,t) < 6, V±.Vg^{v^) ^ 0, wherever Y3Y){x//,x±,t) > 6, 
^ ±-Vg^{v±) < if Ep < -Emax and V ±.Vg^{v±) ~ if = -Emax, and in the small domain 
where Y^d^x//, x±, t) ~ 6, V ±.Vg^{v±) > 0. Hence, in the strongly nonlinear regime the effect 
of V±.Vg^ is mostly dissipative, while it is mostly anti-dissipative in the weakly nonlinear 
regime. However, whatever the nonlinear regime, weak or strong, the anti-dissipative effect 
of ^±-Vg^ always shows up because there is always a region where Fsd > 6 for exiting elec- 
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trons while Ysd < 6 for entering electrons. There, Vgj^{v±) is much larger in magnitude than 
Vg^{—v±) if v±.'V±Ep < 0, so that the effect of Vg^^ is defocussing. In the weakly nonlinear 
regime, defocussing is effective near the center of the wave packet, where the EPW ampli- 
tude is largest. However, as the wave grows, it is clear that the defocussing region moves 
towards the edge of the wave packet and, therefore, becomes less effective. 

It is noteworthy that an analysis similar to ours, on the net energy taken by the electrons 



from the wave due to trapping, was already made in Section IV of Ref. [23| by Banks et 
ai, who insisted on the fact that coUisionless dissipation was in competition with the EPW 
focussing due to wave front bowing. Namely, they wrote that "the focussing effect of curved 
wave fronts is limited by (...) the loss of energy field to resonant electrons that transit 



the wave. . . " . However, in Ref. 



23|, coUisionless dissipation was identified with Landau 



damping while, here, we show that in the weakly nonlinear regime, coUisionless dissipation 
may entai 



a defocussing effect for the EPW, independently of Landau damping. Note that, 
in Ref. [23| , was also discussed the linear defocussing of the EPW due a diffraction-like effect 
which we do not address here. This issue is postponed to a forthcoming paper, where we 
derive an envelope equation for the EPW at second order in its space variations. 

B. Results from simulations of stimulated Raman scattering using the envelope 
code BRAMA 

In order to quantify the defocussing effect due to coUisionless dissipation and to compare 
it with the focussing induced by wave front bowing, we now present results obtained with 
the envelope code BRAMA for an EPW driven by the optical mixing of two counterpropa- 



gating lasers. The code BRAMA, described in detail in Ref. 18[, was recently generalized 



to allow for a three-dimensional space variation of the wave amplitudes, and, as regards 
the electromagnetic fields, it just solves the equations of Ref. [3] with an additional term 
accounting for diffraction. As for the electrostatic wave, its envelope equation is just a mere 
generalization of Eq. iQ obtained by allowing for the coupling to the electromagnetic waves. 
Then, the equation for the EPW amplitude is 



E. 



dtE, + VgVE, + uE, = — ^, (14) 
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where Ed accounts for the effect of the laser drive and where we denoted 

d^xll = j f{v±)d^xfdv^. (15) 

We restrict here to paraxial wave propagation along the x-direction, and approximate 
V±Ep by dyEpj) + d-^EpZ. Then, Vg.VEp = Vg^d^Ep + Vg^dyEp + Vg^d^Ep with 

" "* " K^S '''' 



'"ay,. 



kdu^Xm 



where the transverse components of the wave vector are derived using the consistency relation 
(see for example Ref. 



d,k = -Vu = -^VEp. (18) 



The nonlinear dependence of u with respect to Ep may be found in Refs. 281, |29| , showing 
that du/dEp is usually negative. Hence, the first term in the right-hand side of Eq. ( |T7I) is 
directed along the transverse gradient of the EPW amplitude, i.e., towards the inside the 
wave packet. Therefore, the nonlinear frequency shift, which entails the wave front bowing 
(see Ref. [3l), tends to let the plasma wave packet self-focus. 

In this paper, we restrict to an envelope equation at first order in the space derivatives of 
the EPW amplitude. The inclusion of a diffraction-like effect, as that discussed by Banks et 
al. in Ref. j23|, is postponed to a forthcoming paper where we will allow for nonlinear terms 
proportional to the second order space derivatives of the EPW amplitude. These terms 
are however not essential for a plasma wave driven by stimulated Raman scattering (SRS), 
considered in this Subsection, because an electrostatic wave is much less subjected to diffrac- 
tion than electromagnetic ones since it has a much smaller group velocity. Discussing this 
in detail is however way beyond the scope of this paper, and will be addressed elsewhere. 

Figs. E] to [6] show results from a two-dimensional simulation of the optical mixing between 
a laser with intensity = 6 x lO^^W/cm^ and a counterpropagating "seed" with intensity 
/s = 6 X lO^W/cm^, inside a plasma whose electron temperature is Tg = 0.7keV and electron 
density is = 4 x lO'^ric, where is the critical density. The simulation domain is 140/im 
long and 14/im wide along the ?/-direction. The laser wavelength is 0.527/im and, at best 
focus (located in the middle of the simulation box), its envelope is a Gaussian with intensity 
varying as exp(— y^/w^), with a waist w = Ifim. The seed wavelength is = 0.689/im and 
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FIG. 3: (Color online). Results obtained with BRAMA, at time t = Ips, showing the profiles of: 
panel (a), <1> = eEp/kT^. ; panel (b), the transverse component of the group velocity, Vg^, induced 
by collisionless dissipation, and normalized to the modulus of the group velocity, Vg ; panel (c), 
the total transverse component, Vg^, of the EPW group velocity, normalized to Vg ; panel (d), 
the collisionless damping rate normalized to its linear value. At time t = Ips, when the Landau 
damping rate is still significant compared to its linear value, Vg^ ^ Vg^ , and the sign of Vg^ leads to 
an effective defocussing effect. 



its phase is uniform at the entrance of the simulation box, which is discretized using 1000 
points in the x-direction and 128 points in the y-direction. The distribution of the transverse 
electron velocity is assumed to be a Maxwellian, which we discretize over 25 evenly spaced 
values from —Avth to 4:Vth, where Vth is the thermal velocity. 

As shown in Figs.[3]and|ll in the weakly nonlinear regime, when the Landau damping rate 
is still significant compared to its linear value, the transverse component of the group velocity 
induced by collisionless dissipation, i.e. the last term in the left-hand side of Eq. f ll7p which 
we denote by Vg^, entails a defocussing effect, in agreement with the discussion of the previous 
Paragraph. Moreover, it is easily checked that Vg^ is about ten times larger than the linear 



speed of defocussing due to diffraction-like effects, which, from the formulas of ReL |23| . 



would be t'i™ ^ 8 X 10 ^Vg at time t = Ips. We also note that defocussing is effective near 

^ since ii^max is the same 



'ay 



the center of the wave packet, although, exactly at the center, v] 



9y 
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FIG. 4: (Color online). Same as Fig. [3] but at time t = 1.5ps. As in Fig. [3l the defocussing effect 
due to collisionless dissipation is dominant. Note that, in spite of this, the transverse profile of the 
plasma wave packet has shrunk compared to the previous Figure because of the inhomogeneity of 
the SRS growth rate. 
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FIG. 5: (Color online). Same as Fig. [3] but at time t = 2ps, when Landau damping has nearly 
vanished. At this time, the focussing effect due to wave front bowing prevails near the axis and at 
the front side of the wave packet. 
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FIG. 6: (Color online). Same as Fig. [3] but at time t = 2.3ps, when the component Vg^ due to 
dissipation has changed sign close to the axis and at the front side of the wave packet, which keeps 
on focussing more rapidly. Note the change in shape of the EPW profile between Figs. [5] and [6l 
resulting from the focussing of the front side. 

for entering and exiting electrons. In Figs. E] and SJ the total transverse group velocity, Vg^ 
has the same orientation as Vg^, showing that the defocussing effect induced by Vg^ prevails 
over the focussing effect due to wave front bowing. In spite of this, the transverse extent 
of the plasma wave packet decreases because the SRS growth rate is larger close to the 
axis than near the edge of the simulation box. As the wave keeps on growing, the regime 
of wave-particle interaction becomes more nonlinear and the Landau damping rate keeps 
on decreasing where the EPW amplitude is significant. Then, as discussed in the previous 
Paragraph, the defocussing region moves towards the edge of the plasma pulse so that, 
eventually, the focussing effect due to wave front bowing dominates near the center and at 
the front side of the plasma wave packet, as shown in Fig. [51 This clearly affects the profile 
of the EPW, as may be seen by comparing the results of Fig. [5] to those of Fig. [61 Finally, 
in the strongly nonlinear regime, when u ^ where the EPW amplitude is significant, Vg^ 
changes sign close to the axis and on the front side of the wave packet, which therefore keeps 
on strongly focussing. 



19 



IV. CONCLUSION 



In this paper, we discussed in detail the distinction that may exist between Landau 
damping and colhsionless dissipation. Three regimes of wave-particle interaction actually 
need to be considered. When the electron response to an EPW is well approximated by 
linear theory, the wave is Landau damped, its amplitude decreases along the characteristics, 
while its group velocity remains constant and uniform. In the linear regime, dissipation is 
only due to Landau damping. In the nonlinear regime and in one dimension, the EPW group 
velocity, Vg, remains nearly at its linear value as long as F = ojb{x — v^t', t')dt' < 6, and 
quickly increases towards the phase velocity where Y i.e. where the trapped electrons 
have completed about on bounce period, which entails a local increase in the size of the wave 
packet. For larger values of F, Vg decreases with x (i.e. along the direction of propagation of 
the wave packet), until the maximum wave amplitude is reached, and then remains constant. 
This clearly induces the shrinking of the wave packet in that region located between x = xq, 
where y 6, and x = Xmax where the EPW amplitude reaches its maximum value. In 
the weakly nonlinear regime, when Y ^ Q close the where the EPW 

amplitude is largest and then remains constant for x > Xmax- Hence, the nonlinear variations 
of Vg lead to an increase in the size of the wave packet, which is an anti-dissipative effect. As 
for the colhsionless damping rate, it remains nearly at its linear value for x < Xmax so that, 
again, damping is the only dissipative mechanism. By contrast, in the strongly nonlinear 
regime, when y 6 in the rear side of the wave packet where the wave amplitude is very 
small compared to its maximum value, only there is the damping rate significant compared to 
its linear value, so that damping is much less effective than in the linear or weakly nonlinear 
regimes. As for the nonlinear variations of Vg, they lead to the shrinking of the wave packet 
in the region xq < x < Xmax, where the EPW amplitude is significant, which automatically 
reduces the electrostatic energy W = E^dx. This a dissipative mechanism different from 
Landau damping, and more effective for longer and more intense wave packets. Physically, 
it is directly due to trapping, which is dissipative, since the energy gained by the electrons 
while being trapped is not given back to the wave when they are symmetrically detrapped. 
In Ref. (2l|, we wrote that Vg ^ dtu/dk due to a term in the envelope equation that may 
be viewed as the nonlinear counterpart of the Landau damping rate, although the result 
Vg 7^ du/dk was only valid once Landau damping has nearly vanished. Hopefully, in this 
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paper, we clarified this statement by showing that the nonhnear variations of the EPW group 
velocity, as derived in Ref. j21], was one manifestation of colhsionless dissipation effective 
in the strongly nonlinear regime, i.e., after Landau damping has nearly vanished where the 
wave amplitude is significant. 

The other important point of this paper was to recognize how transparent was, in the 
envelope equation for the EPW amplitude, the effect of trapping on colhsionless dissipation. 
Once this is understood, the envelope equation derived when the wave amplitude is uniform 
is straightforwardly generalized to allow for three-dimensional effects. The conclusions are 
essentially the same as in one dimension. In the weakly nonlinear regime, dissipation is 
only due to Landau damping, while the increase of the group velocity compared to its linear 
value would induce an anti-dissipative effect by increasing the space extent of the wave 
packet. In the strongly nonlinear regime, damping is strongly reduced, and the nonlinear 
variations of the group velocity entails a shrinking of the wave packet and, therefore, of 
the electrostatic energy. However, we came to the quite unexpected result that dissipation 
induced a transverse component to the EPW group velocity, which we precisely quantified. 
In the weakly nonlinear regime, this transverse component is directed towards outside the 
wave packet, leading to a defocussing effect that might overcome the focussing due to wave 
front bowing, as illustrated by results from two-dimensional simulations of optical mixing 
using the envelope code BRAMA. As the EPW amplitude keeps increasing, and the wave- 
particle interaction becomes more non linear, the transverse component of the group velocity 
due to dissipation eventually entails a focussing effect that reinforces that due to bowing. 



At this stage, we would like to point out that, as clearly shown in Refs. 2^, |25|, symmetric 
detrapping leads to a positive slope in the electron distribution function, that may entail 
the unstable growth of electrostatic modes. This may be viewed as an indirect consequence 
of colhsionless dissipation, which will however not be discussed in this paper. This is left 
for future work. 
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